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Abstract 

Exactly solvable two-dimensional polygon models, counted by perimeter and area, 
are described by g-algebraic functional equations. We provide techniques to extract 
the scaling behaviour of these models up to arbitrary order and apply them to some 
examples. These are then used to analyze the unsolved model of self-avoiding poly- 
gons, where we numerically confirm predictions about its scaling function and its 
first two corrections to scaling. 
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1 Introduction 

Two-dimensional polygon models [521 (T5| such as self-avoiding polygons f^lEllini, counted 
by perimeter and area, are described in the grand-cononical ensemble where one introduces 
a perimeter activity x and an area activity q. The perimeter and area generating function 
G{x, q) is defined as 

oo 

G(x,g) = Y,Vm,nx'^q^ = 5^p™(g)x'", (1.1) 

m,n ni=0 

where Pra,n denotes the number of polygons on a given lattice with perimeter m and area 
n. The thermodynamical properties of polygons are described by the behaviour of the 
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perimeter and area generating function about its singular points. The singular behaviour 
of the generating function translates to properties of the coefficients Pm,n in the limit of 
large perimeter or area. The locus of singularities (phase diagram) of a typical polygon 
model is shown in Fig. 1. 
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Figure 1: Phase diagram of a typical two-dimensional polygon model 



Of particular interest are points where lines of singularities meet, which describe phase 
transitions. For example, the model of self-avoiding polygons (SAPs) is believed to exhibit 
tricritical behaviour, to be described below, about the point where the perimeter generating 
function 1) is singular. This point describes a crossover between extended and deflated 
self-avoiding polygons: The mean area (a)m of SAPs of perimeter m grows asymptotically 
like [131 

\A(g)m (0<g<l). ^^'^^ 

It can be shown J3] that in the limit g — >^ the generating function is dominated by 
polygons of minimal area. Since these polygons may be viewed as branched polymers J2] , 
the phase g < 1 is also referred to as the branched polymer phase. 

Little exact information is known about SAPs El 13 El? almost all information being 
as a result of numerical studies p[Hl UH UH CE]- Recently, a prediction of the functional 
form of the scaling function, which describes the tricritical behaviour, has been derived and 
conflrmed numerically |29j . This led, by re-interpretation using fleld-theoretic methods, to 
the prediction of a number of other exact scaling functions [S]. In this paper, we explain 
our previous analysis j2H] within a larger mathematical framework and extend it to include 
analysis of corrections to scaling. 

Due to the difficulty of the SAP problem, simpler polygon models, which are subclasses 
of self-avoiding polygons, have been analyzed and solved exactly [23 121 ESI- Over the years, 
it has become clear that the underlying structure is described by g-algebraic functional 
equations [H UHl HI ISl 121 HHj , which we will analyze in detail in this paper. Assume that the 
perimeter and area generating function G(x, q) of a polygon model satisfies a g-algebraic 
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functional equation [21] 

P{G{x, q),G{qx, q),..., ^(g^x, q),x, q) = 0, (1.3) 

where P(?/o; Uii ■ ■ ■ i Vn, q) is a polynomial in ?/o; Vi, ■ ■ ■ , Vn, x and q. The motivation for 
this type of equation arises from the requirement that polygons of a given size can be built 
up recursively from polygons of smaller size [HI 123 l2| • The existence of such a recursion 
accounts for the solvability of the model. We call polygon models of such a form g-algebraic. 

The limit g ^ 1 in p.3|l leads to an algebraic differential equation for the perimeter 
generating function G{x, 1) of order A^, i. e. a vanishing polynomial in x, G{x, 1) and 
its derivatives. From g-linear equations, the limit q —>■ 1 leads to a linear differential 
equation of order N. This may be seen by first expressing G{qx, q) in terms of the q- 
derivative {DqG){x) = {G{x, q) — G{qx, q))/{l — q)x, likewise G{q^x, q) in terms of higher 
g-derivatives, and then expanding the functional equation about g = 1. The order (1 — g)" 
in the expansion contains derivatives of degree not exceeding n. The defining equation for 
the perimeter generating function G(x, 1) is given by the non-trivial contribution of lowest 
order. One simple case is given by n = 0, where derivatives are absent. In this case, the 
perimeter generating function satisfies the algebraic equation 

P(G(x, 1), 1), . . . , G{x, 1), X, 1) = 0. (1.4) 

Generally, exactly solved polygon models have algebraic perimeter generating functions 
[2]. However, it has been shown [27] that the anisotropic perimeter generating function 
of the unsolved model of self-avoiding polygons is not algebraic (even not D-finite |30j). 
and there are other simple, unsolved polygon models such as three-choice polygons [TT] . 
which displays a logarithmic singularity. A natural question to ask is if these models are 
also described by g-algebraic functional equations. As we showed ^2^, for rooted self- 
avoiding polygons G^^\x,q) = x-^G{x,q) this assumption appears to lead to the correct 
prediction for the singular part of the perimeter and area generating function about the 
phase transition point, thereby answering a conjecture by Guttmann about the scaling 
function ^23] in the affirmative. Of course, this does not imply that rooted self-avoiding 
polygons are g-algebraic. 

From physical grounds, we expect tricritical scaling behaviour about such points, 

G(x, g) ~ G^^'^\x, g) + (1 - qYF{{x, - x){l - g)"*), (x, g) (x,", 1"), (1.5) 

where F{s) is a scaling function, depending only on a combined variable s = (xc — x)/(l — 
g)'^, which is commonly assumed to be regular at the origin. Its value at the origin deter- 
mines the singular behaviour of the area generating function G'^**"^)(xc, g) ~ F(0)(1 — g)^ 
about g = 1. The behaviour of the scaling function at infinity determines the singu- 
lar behaviour of the perimeter generating function G{x, 1) about the critical point Xc- If 
F{s) ~ as^'^ for s ^ oo, it follows that G{x, 1) ~ a(xc — x)"'^ for x Xc, if 6' + 07 = 0. 
The exponent (p is called the crossover exponent. The slope of the critical line Xc(g) about 
g = 1 is determined by the first singularity Sc of F{s) on the negative real axis. We have 
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Xc{q) ~ Xc — Sc(l — q)'^ about q = I- If F{s) ~ a(s — Sc)° about s = Sc, the singular part of 
the perimeter and area generating function behaves as G^'^^'^^\x, q) ~ a{l—q)~°"^{xc{q)—x)°' 
for X — >• Xc{q). 

From the scahng function of the grand-canonical ensemble, finite-size quantities may 
be computed. It can be shown that for large m under mild assumptions 



Prn{q) ~ m'-^x-r-'h ( -(1 - q)-^ ] (m ^ oo, g ^ 1"), (1.6) 



where the finite-size scaling function h{x) is given by 

oo „ 

(1-7) 

and the coefficients fk appear in the asymptotic expansion of the scaling function 

oo „ 

^W = Esra)«- (1-8) 

fe=0 

Due to the appearance of the Gamma function in ()1.7|1 . the functional form of the finite-size 
scaling function is generally more complicated than the grand-canonical scaling function. 
The determination of scaling functions using an exact solution is known to be a difficult 
problem. They have been extracted for semicontinuous versions of a number of polygon 
models such as Ferrers diagrams and stacks [25, staircase polygons [231122] and column- 
convex polygons [Hj. In some cases, explicit expressions for the finite-size scaling function 
can be derived |2I|. 

Recently it has been shown that the scaling function for simple discrete models can be 
derived directly from the defining g-functional equation [211122] 5 and one is led to ask about 
the generality of the approach. This is the aim of the present paper. As a first step towards 
understanding the scaling behaviour of general g-algebraic polygon models, we focus on 
the case where the perimeter generating function is algebraic. Remarkably, the generic 
form of the scaling behaviour of the perimeter and area generating function is determined 
by the nature of the algebraic singularity alone, irrelevant of the details of the functional 
equation. This we illustrate in the simplest case of a square-root singularity and a simple 
pole at X = Xc- In these cases, we will perform a formal asymptotic expansion about 
X = Xc and g = 1 which is uniform in q. This asymptotic expansion does not only give the 
scaling behaviour but also gives corrections to the scaling behaviour up to arbitrary order. 
After presenting exactly solvable examples, we discuss the model of self-avoiding polygons. 
The perimeter generating function of rooted self-avoiding polygons displays a square-root 
singularity, but is almost certainly not algebraic [27] . We will confirm numerically that 
the scaling function about the singular point is the one predicted by the assumption of a 
g-algebraic functional equation for the perimeter and area generating function. We will 
also confirm predictions about the first two corrections to scaling. The paper concludes 
with a discussion of open questions and possible future work. 
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2 Area moments of (/-algebraic polygon models 



Let G{x, q) be the perimeter and area generating function of a polygon model. The area 
weighted moments gk{x) appear in the formal expansion about q = 1 

G{x,q) = Y,Pm,nX^q- = - q)' (2.1) 

m,n k 

and are of the form 

^'^^^ " "^Tr 5^ ^(^ - 1) ■ ■ ■ - ^ + (2.2) 

m,n 

Note that all area moments have the same radius of convergence as the perimeter generating 
function go{x). This can be seen by a ratio test, using the estimate 

n{n — 1) ■ ■ ■ {n — k + 1) < n'' < m^'' (n>no). (2.3) 

The second inequality amounts to the fact that squared perimeter is larger than the 
area n of a polygon. 

Let us assume that a polygon model satisfies the g-algebraic functional equation 

Piyo,yu---,yN,x,q) = 0, yj, = G{q''x,q). (2.4) 

We will show below that the functional equation determines the moments gk{x) recursively. 
In order to proceed, we have to specify the properties of the functional equation in the 
limit g — i> 1. Therefore assume that the limit q ^ 1 leads to the algebraic equation 

^'(2/0,1/0, • • • ,Z/o,a;, 1) = 0, yo = G{x,l). (2.5) 

The singular behaviour of the perimeter generating function G{x, 1) is determined by the 
algebraic equation ()2.5p [T^ . This in turn determines the singular behaviour of all other 
moments. We will make this explicit for the two simplest forms of algebraic singularity: a 
square-root singularity and a simple pole. This generalizes previous investigations of the 
model of staircase polygons, where singular exponents of the first two area moments have 
been determined 

2.1 Recursive determination of the area moments 

We show that the area weighted moments can be computed recursively from the functional 
equation if the limit g — 1 leads to a non-trivial algebraic equation for the perimeter 
generating function go{x). 

It is instructive to first derive an expression for gi{x). To this end, we expand the 
g- algebraic functional equation to first order in e = 1 — g, making use of the expansion 

Giq^x, q) = goix) + eigiix) - kxg'^{x)) + 0{^). (2.6) 
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The terms of order e are given by 

In the above equation, P is evaluated at argument (?/o, . . . , Un^ x, q) = {yo, . . . ,yo,x,l), 
and we adopted the notation dk = dy^ . The first term is non-zero by the assumption of an 
algebraic equation for go{x) of the form ()2.5|1 . This means that gi{x) is given in terms of 
powers of gQ{x) and its derivative. In the general situation the following statement holds: 

Lemma 1. The order e" in the expansion of the functional equation is linear in gn{x) and 
contains products of gk{x) and its derivatives with Q < k < n. 

Proof. This is seen by an analysis of the terms appearing in the expansion of the functional 
equation in powers of e. We will do this in detail for subsequent purposes. Note first that 
the expansion of a g-shifted function f{q^x) is given by 

fiq'x) = J2i-^y ( ^ ' / W = E E ^' ^)^"^/^'"^ (^), (2-8) 

1=0 ^ ^ Z=0 m=0 



h{k, l,m)=^-^[J2 Si-JS^lF ] . (2.9) 



where 

Here, Sm,n and Sm,n denote the Stirling numbers of first and second kind, respectively. We 
have h{k, 0, 0) = 1 and h{k, /, m) = for m > /. The expansion of G{q^x., q) is then given 

by 

G{q^x,q)=gQ{x) + hk{e,x) (fc = 0, . . . , iV). (2.10) 
The function hk{e,x) has an expansion in powers of e of the form 

oo n I 

h{e, x) = J2 ^''hk,n{x), hu,n{x) = E E ^' ^)^"&{^)- (2.11) 

n=l l=Q r=0 

With the introduction of 

i7= (/io(e,x),...,/;,Ar(e,x),-e), V = (9o, . . . , Sjv, 9g), (2.12) 
the expansion of the g-functional equation can be written as 



oo 



P{yo, ■■■,yN,x,q) = ^^^P(yo, • • • , 2/0, a;, 1) 

fc=0 

oo (— e)'' 



r,m=0 k„ 



(2.13) 



Note that the first sum is finite since P{yo, • • • , Vn, x, q) is assumed to be a polynomial in 
its variables. In the second equation we introduced the abbreviation x„ = xi, . . . , x„. The 
range of summation is < fcj < A^. We now insert p.llj) into the above equation and 
extract the term of order e". It is given by 

E E E-""' ff[M*..«..*.)9;:i.(:^)) • 

r,m=0 ■ ■ \rm\ ^™ \i=l / (2-14) 

=n—r ^™ 

■<9g<9k„P(?/o, •• -,1/0,3;, 1), 

where we introduced the abbreviation |x„| = xi + ■ ■ ■ + x„. We have 1 < Tj, < Sj < rj and 
< tj < Sj. The summations over and are independent. It follows from ()2.14|) that 
the order contains gn{x) only linearly: The constraint Ir^l = — r implies r = and 
m = 1. This in turn implies Si = ti = and leads to gn{,x) dkP{yo, . . . ,yo,x,l). As we 
argued above, this term is non-zero. The other terms are given by products of derivatives 
of gk{x) with k < n. It is therefore possible to solve for the area moments recursively. □ 



2.2 Behaviour about a square-root singularity 

Assume that the algebraic perimeter generating function G{x, 1) has a square-root singu- 
larity at X = Xc. It follows from the theory of algebraic functions ^1] that this is equivalent 
to the conditions 

(Y. -^fc^) = 0' {Y1 ^fc.'^) ^ 0' (^--P) ^ 0' (2.15) 

on the functional equation, evaluated at x = Xc. We have introduced the abbreviations 
9k = dy^ and dk^i = dyf^dyi for < k,l < N. The first condition implies that x = Xc is 
a singular point, the second and third condition indicate a square- root singularity. The 
expansion is of the form 

oo 

go{x) = goixc) + ^ fo,iixc -x)~. (2.16) 

1=0 

The coefficients /o,; can be determined recursively from the algebraic equation (|1.4|) . using 
the constraints ()2.15|) . We have in particular 

flo = . ^^^"^^ V (2.17) 
l{Ek,idk,iP) 

The sign of /o,o determines the branch of the solution. Since we are dealing with generating 
functions having non-negative Taylor coefficients, the negative sign has to be chosen. With 
the expressions given above it is easy to analyze the singular behaviour of gi{x). To this 
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end, we differentiate the algebraic equation ()1.4|) in order to replace (^j. dkP) in ()2.7|) . We 

get 



We now make the assumption that 



J2kdkP]^0. (2.19) 

V k / 



This generic case is satisfied in the models discussed below. It then follows that gi{x) 
admits a series expansion of the form 



oo 



^7i(a;) = Xl/i.'«(^'=-^)^"'' (2-20) 

fc=0 

where the coefficients fi^k can be computed from ()2.18p . The leading coefficient /i,o is 
explicitly given by 

(y,.kdkP) 

-4/1,0 = X, y^'^ \ . (2.21) 
-2[E,,id,,iP) 

For moments of higher order, we have the following result: 
Theorem 1. 

Under the assumptions (j2.15p and (|2.19p . the area moment gn{x) admits an expansion 
about its singular point Xc of the form 

00 

9n{x) = fnAxc - x)^-i"+i (2.22) 

k=0 

It can be computed recursively from the g-algebraic functional equation at order e". The 
coefficients /„,o are explicitly given by 

fnfl = Cn/i, 0/0,0 ' (2.23) 

where the coefficients /i,o and /o,o are given by ()2.21|) and ()2.17|) . and Cn satisfies the 
quadratic recursion relation 

^ n— 1 

C„ + (3n - 4)c„_i + ^ Cn-rCr, Cq = 1. (2.24) 

Remarks: 

(1) The behaviour of the leading coefficients is independent of the degree M > 2 of the 
g-algebraic functional equation. 

(2) The first few values of the coefficients are Ci = 1, C2 = —5/2, C3 = 15, C4 = —1105/8, 
C5 = 1695, Cq = -414125/16, Cj = 472200, Cg = -1282031525/128, Cg = 242183775, 
Cio = -1683480621875/256. 
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Proof. We first prove the stated form of the expansion of the functions gn{x) about the 
singular point. This we do by induction. Since the case n = 1 has been treated above, 
assume now that the statement holds for gk{x) where k < n. We analyze the g-algebraic 
functional equation at order e". By Lemma 1, it contains gn{x) linearly, and the other 
terms contain products of derivatives of gk{x) with k < n. From these facts it follows 
that X = Xc is a singular point of gn{x) and that there exists an expansion of gn{x) about 
X = Xc in half-integer exponents. To get the leading singular exponent of gn{x), we analyze 
the singular behaviour in a typical term of the form ()2.14|) . which is determined by the 
product appearing there. Terms not containing gn{x) have by induction assumption a 
negative exponent of modulus 

^^n-r)-- + Y,[U-^sA. (2.25) 

i=l ^ ' 

The leading singular behaviour is obtained by its maximum value subject to the constraints 
< ti < Si. The maximum value appearing is 3?2/2 — 1. It is realized by g'^_i{x) with 
values (r, m, ti, si) = (0, 1, 1, 1) or by terms quadratic in gk{x) with values {r,m,ti, Si) = 
(0,2,0,0). This leads to 

= 0((x, -x)i-i"). 

We can read off the leading singular behaviour of gn{x) from this equation. Since the prefac- 
tors appearing are all non- vanishing, we get an exponent of (1 — 3?t,)/2. The above equation 
also determines the values of the leading coefficients ()2.24j) . as is seen by a straightforward 
computation. □ 

2.3 Behaviour about a simple pole 

In this case it is advantageous to rewrite the g-algebraic functional equation in the form 

M 

P{yo,---,yN,x,q) = ^^ak„(x,g)t/fci---?/fc„ + ao(x,g). (2.27) 

m=l k„ 

We introduce the abbreviation am{x,q) = '^k™(a;, g) and am{x) = am{x,l). In the 
limit g — >■ 1, the perimeter generating function go{x) satisfies the algebraic equation 

aMix)goix)'^^ + aM-i{x)go{x)^'~'^ H h ai{x)goix) + ao(x) = 0. (2.28) 

We assume that at the critical point x = Xc the algebraic equation satisfies 

aM{xc) = 0, a'^ixc) 7^ 0, aM^i{xc) 7^ 0. (2.29) 
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This is equivalent to the assumption of a simple pole. From the theory of algebraic functions 
[T^ it follows that go{x) has an expansion of the form 

oo 

9o{x) = Y,foA^c-x)'-\ (2.30) 

fc=0 

where the coefficients /o,fc can be calculated from the algebraic equation. In particular, the 
first coefficient is given by 

/o,o = ^(Xe). (2.31) 
An expansion of the g-algebraic functional equation to first order in e gives 

M \ /A/ 



9ii^) I ^mam{x)go{x)'^ M ~ I 5Z (^9'^™) (2^, l)5'o(a;)"' 

m=l / \m=0 

g'oi^) U E ( E Ik™|ak„(x, 1) ) go{xrA = 0. 



(2.32) 



m=l \ k„ 



This is exactly ()2.7|) . but the present notation is better suited to analyze the singular 
behaviour of gi{x), since the divergent behaviour of the perimeter generating function 
go{x) has to be taken into account. The above equation shows that the first area moment 
has an expansion about x = Xc in integer powers. We assume that 

(ElkAflak,,) (xe,l)^0. (2.33) 

This generic behaviour is satisfied for the model discussed below. The leading singular 
term in the equation is then of order {xc — a;)"^^~^, and we conclude that gi{x) has an 
expansion of the form 

00 

gi{x) = J2fiAxc-x)'-\ (2.34) 

k=0 

The expansion coefficients can be computed from the above equation. The first coefficient 
is given explicitly by 

r ^ (EkM \^M\akM) QAf-l , . 

—Jl,0 — Xc } — . [Z.6b) 

All expressions are evaluated at argument (x, q) = {xc, 1). In the general case, the following 
statement holds: 

Theorem 2. Under the assumptions ()2.29|1 and ()2.33|1 . the area moment gn{x) admits an 
expansion about its singular point Xc of the form 

00 

9n{x)=Y,fnAXc-x)''-'^^-\ (2.36) 
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It can be computed recursively from the g-algebraic functional equation at order e". The 
coefficients /„_o are given by 

fn,0 = Cn/l,o/o,0 i (2.37) 

where /o,o and fi^ are given by (j2.3H) and (j2.35|) . respectively. The numbers c„ satisfy the 
linear recursion relation 

Cn + {2n - l)cn-i = 0, Co = 1. (2.38) 

Remarks: 

(1) The result is independent of the degree M > of the functional equation. 

(2) The coefficients are given explicitly by c„ = i—'^)^ 2"'^^{n-iy. • "^^^ ^^^^ values are 
ci = -1, C2 = 3, C3 = -15, C4 = 105, C5 = -945, Cq = 10395^7 = -135135, Cg = 2027025. 

Proof. This statement is proved similarly to the square-root case treated above. Assume 
that the statement holds for gk{x) where k < n. We analyze the singular behaviour of the 
contribution from order e". As follows from Lemma 1, the point x = Xc'is a. singular point 
of gn{x) and allows for an expansion of gn{x) in integer exponents. Now consider terms 
containing g'^_^{x). The prefactor appearing is exactly the second bracket in ()2.32p . Using 
the induction assumption, this results in a divergence with exponent 2?2 + M — 1. This is 
also the leading divergence: Consider the singular behaviour in a general term in the sum 
(|2.14p . Using the induction assumption again, we see that the product appearing there 
diverges with an exponent of absolute value 

m 

2 (n - r) + m + ^ (ti - 2si) . (2.39) 

i=l 

The derivatives of P in ()2.14j) contain diverging terms as well. It can be computed that 
the most singular contribution is go{x)^^~"^. If Si = ti = 0, this term is accompanied by 
om, adding an exponent M — m — 1 to ()2.39|) . due to the pole condition. Otherwise, the 
divergence adds an exponent M — m. This leads, by power counting, to two different 
choices of exponents: (r, m, tj, Sj) = (0, m, 0, 0) or (0, m, 1, 1) for exactly one i. An explicit 
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calculation of the prefactors appearing yields the functional equation to leading order 

M 



m=l \|r™,|=n 



ml \ m 



M 



(2.40) 



X 



m=l 



I 9'r,^li^)9r2{x)---gr 



m — 1 



An inspection of the m = 1 term reveals a divergence with exponent 2n + 1 for gn{x), as 
claimed above. From the above equation a recursion relation for the values of the leading 
coefficients can be obtained. It is 



M 

E 



M -1 
m — 1 



|rm,|=n 



+ 



(2ri - l)cri-iCr2 ■ ■ - Cr 

\rm\=n 



(2.41) 



0. 



This equation determines the values c„ recursively, which may be seen by induction on n, 
using > 0. We get the /inear recursion ()2.38p . being independent on M. □ 



2.4 Area moments and scaling behaviour 

We now explain how the singular behaviour of the area moments is related to the scaling 
function and the corrections to scaling. Assume that the expansion of the area moments 
gk{x) about the singular point x = Xc is given by 



kl 



^ (Xc - Xpi^'' 



(2.42) 



where 7fc,;+i < '~ik,h and fk,i are the amplitudes. We next assume that the exponents 7^,/ 
are of the special form 

lk,i = ---■ (2-43) 

This implies that 6*^+1 > 9i and means that we can order the singular terms columnwise, 
where in each column the exponent increases by a constant indexed by k. In this special 
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situation, which is satisfied in our examples, the singular part of the perimeter and area 
generating function has a formal expansion of the form 

oo 

= ^e^'F; {{xc - x)e~''') , (2.44) 
/=o 

where Fq{s) is the scaling function and Fi{s) for / > 1 are the correction-to-scaling 

functions. We note scaling exponents 6q and (f) may be obtained from the analysis of 
the leading singular behaviour of the first two area moments. This technique was used 
intensively before [nil2Sl- The amplitudes fk^i appear in the asymptotic expansion of the 
functions Fi{s) via 

oo J, 
fc=0 

In the next section we will derive differential equations for the scaling and correction-to- 
scaling functions Fi{s) directly from the functional equation. 



3 Scaling of g^-algebraic functional equations 

Assume that the perimeter generating function has an algebraic singularity at x = Xc- As- 
sume further that the singular part of the perimeter and area generating function displays 
scaling behaviour about x = Xc and g = 1 of the form 

Gix, q) ~ G(^-^)(a;, g) + (1 - q^F, {{x, - a;)(l - q)-^) (x, q) ^ (x' , l"). (3.1) 

We expand the functional equation in terms of the scaling function Fo{s) at argument 
s = {xc — x){l — q)~'^. In the limit g ^ 1, this will lead to a differential equation for the 
scaling function, as has been demonstrated for a number of simple examples previously 
[2H1 121] ■ The technique is inspired by the method of dominant balance, which has been 
used to derive scaling functions for semicontinuous versions of polygon models [201 IHl 12^] • 
Introducing e = 1 — g, we first give the expansion of a g-shifted function in analogy to ()2.8|) 

F ((X. - q'x)e-^) = f^^-eY (^^^ ' F(.) 

oo m y \ oo V • / 

= ^i7^M(s)^ (-a;e)"("')s'"""5^e'-"*/i(A;,/,m). 

m=0 n=0 ^ ^ l=m 

This implies that the m-th derivative scales with an exponent I — n(j) in e, where < 
n < m < I. The smallest exponent occurs for I = n = m. In the expansion of the scaling 
function of g-shifted argument this leads to m-th derivatives F^\s) with scaling exponent 
9 + m(l — 0). If > 1, this restricts the scaling function to be a polynomial in s: In the 
expansion of the functional equation, derivatives of the scaling function will appear with 
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arbitrarily small exponents in e. This means that the limit e ^ can only exist if almost 
all derivatives of the scaling function vanish, i.e. if the scaling function is a polynomial. 
We assume in the following that < 1. This assumption is satisfied for the cases treated 
below and will result in a differential equation of first order for the scaling function. If 
0=1, higher derivatives cannot be ignored, leading to a difference equation for the scaling 
function [2H]- The expansion of the generating function with g-shifted argument is given 
to lowest order by 

^(g^^x, q) = G^'-'s^x, q) + F^{s) + kx.e'^^^-'*'^ F^{s) + 0(e^+2{i-<^))_ ^3 3) 

We now analyze the case of a square-root singularity aX x = in detail. 



3.1 Scaling function about a square-root singularity 

We first derive a differential equation for the scaling function about a square-root sin- 
gularity and solve it. In order to keep notation simple, we omit arguments, which are 
understood to be always at (7/0, . . . , Un, x, q) = {yo, . . . ,yo, Xc, 1). Using ()3.3j) . we expand 
the functional equation ()1.3|) in powers of e. We assume a square-root singularity, which 
implies (p = 26 and 



\ fc / \ k,i 



P = 0, dkPj =0, ( > ^ dk,iP I ^ 0, {d,P) ^ 0. (3.4) 

The lowest orders in e are given by 



5^fcPjF^(.) + --- 

k ) 



(3.5) 



The assumption that the two exponents do not coincide leads to a trivial scaling function. 
In order to have a nontrivial scaling function, we therefore demand 2^ = = 6'+(l — 0) 
and get exponents 

Q = \ 0=1 (3.6) 

This coincides with the derivation using the area moments. We obtain a differential equa- 
tion for the scaling function of the form 

Fo(s)'-4/i,oF,;(s)-/o> = 0, (3.7) 

where the coefficients /i,o and /o,o are given by 

2 {d^P) „ _ {Y^^kd^P) ...^ 

/o,o - 7 V' •'I'O ~ ~"V7 
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The solution is uniquely determined by the prescribed asymptotic behaviour at infinity 
and is given by 

f„W = -4/4l.A,((AL)^'%j. (3.9) 

where Ai(x) = ^ cos(t^ /3 + tx) dt is the Airy function. The coefficients fkfl in the 
asymptotic expansion of the scaling function are exactly those derived in the previous 
section. The asymptotic behaviour is given by 

Fo{s) ~ hfls^'^ (s^oo) (3.10) 

f f 

The first singularity on the negative real axis is a simple pole. It is located at ( j Sc = 



3.2 Corrections to scaling about a square- root singularity 

The analysis of the singular behaviour of the area weighted moments justifies a formal 
scaling Ansatz 

oo 

G(-"^)(a;,g) = J]e('+i)^F, ((x, - a;)e-^) . (3.11) 

l=Q 

We show that the order e'-'"'"^''^ term in the expansion of the functional equation yields 
a linear, inhomogeneous differential equation of ffrst order for the correction-to-scaling 
function Fi{s). First note that, due to ()3.2|) . the expansion of the g-algebraic functional 
equation contains terms of the form 

with non- negative integers ri, r2, Vij. Here, f/'^^ denotes the k-ih. derivative of the 
correction-to-scaling function -F)(s), and we have the constraints < -n,-,- < kj < mj. 
Let us analyze the order e^'"*"^^^ in detail. It consists of terms quite similar to the expansion 
in ()3.5p . The contributions from the highest correction-to-scaling functions are 

J29kP] Fi+^{s)+x, ($^A;afcP) F;{s)+ (5^9fc,/P) Fi{s)Fois). (3.13) 

\ k / \ k J \k,l / 

Since the first prefactor vanishes by the assumption of a branch point, F;_|_i(s) does not 
contribute to that order. There are however non- vanishing contributions linear in Fi{s) and 
its derivative. The remaining terms contain products of the scaling function and corrections 
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to scaling of lower order than / together with their derivatives. Given the functions Fk{s) 
for < / and the desired asymptotic behaviour as discussed in the previous paragraph, 
integration of the above equation uniquely determines Fi{s). 

We note that for special solvable cases like staircase polygons, the corrections to scaling 
have a simple form. They can be written as a polynomial in the scaling function and the 
scaling variable. This is demonstrated below. 



3.3 Scaling function about a simple pole 

The simple pole is characterized by 

aM{xc) = 0, a^(xc) 7^ 0, aM-iix^) ^ 0. (3.14) 

About a simple pole, the perimeter generating function is diverging. This means that there 
is no regular term in ()3.1|) and that the exponent 6 is negative. Let us therefore write the 
leading behaviour about the critical point as 

G{x, g) ~ (1 - q)-'Fo{{x, - x)(l - g)-^) (x, q) ^ (x,", 1"). (3.15) 

Note that we have redefined 6 such that it is positive. Since we are dealing with a simple 
pole, we can assume 6 = (p. We now derive a differential equation for the scaling function 
Fq{s). This is done by expanding the functional equation about s = {xc — x)e~'^ and taking 
the limit e — > 0. To this end, we first multiply the functional equation by the exponent of 
its leading divergence e^'^^, where M is the degree of the functional equation. To lowest 
order in e, we get contributions from order M and order M — 1 of the form 



= -a'^,sFo{s)^'e^ + I \^mK,s ) x,F;,{s)Fo{s 



+ aM-iFo(s 



M 

M-1, 



(3.16) 



In order to obtain a non-trivial scaling equation, we demand equality of the three expres- 
sions and get exponents 

= 1, ^=\- (3.17) 

These exponents follow already from the analysis of the area weighted moments of the 
previous section. Here, their values arise from a consistency argument. Assuming that the 
second term in the above expansion does not vanish, we get a differential equation of the 
form 

sFo(s) + ^F^(s)-/o,o = 0, (3.18) 
/o,o 

where the coefficients /i^o and /o,o are given by (j2.31|) and (j2.35|) . The solution of this 
equation is uniquely determined by its prescribed asymptotic behaviour and is given by 



Fo{s) = /o.ov^erfc {^s) e''' , (3.19) 
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where he have introduced the abbreviation a = — /o,o/ (2/i,o), and erfc(x) = e~* dt 

is the complementary error function. The coefficients /„^o in the asymptotic expansion of 
the scahng function are exactly those derived in the previous section. The asymptotic 
behaviour is given by 

Fo{s) ~ fo,os-' (s^oo) (3.20) 

^o(O) = /o,oV™ 

Fo{s) ~ 2Fo(0)e"^' (s ^ -oo) 

The first singularity of the scaling function on the negative real axis is at minus infinity, 
which implies Xc{l) > Xc for the critical line Xc{q) 



3.4 Corrections to scaling about a simple pole 

The analysis of the singular behaviour of the area weighted moments justifies the formal 
scaling Ansatz 

oo 

G{x, q) = Yl ^"^'''^'Fi ((^c - x)e-^) . (3.21) 

1=0 

As in the previous paragraph, we multiply the functional equation by in order to 
obtain a finite limit as e — > 0. We show that the order e^'"*"^-*^ term in the expansion of the 
functional equation yields a linear, inhomogeneous differential equation of ffist order for 
the correction to scaling function Fi{s). This is done in analogy to the case of a square- 
root singularity discussed above. Note ffist that the expansion of the g- algebraic functional 
equation leads to terms of the form 

with non- negative integers ri, r2, Vij. We have the constraints < rij < kj < rrij. Let us 
analyze the order e('+^)^ term in detail. The terms with highest corrections to scaling are 



<-M 

- 2a'MsFi{s)Fois)^'-' 



AI-l 

(3.23) 



The ffist term does not contribute since the prefactor vanishes by the assumption of a 
branch point. The following terms contain Fi{s) and its derivative linearly. The remaining 
terms contain products of the scaling function and correction-to-scaling functions of lower 
order than / together with their derivatives. Given the functions Fk{s) with k < I and the 
desired asymptotic behaviour as discussed in the previous paragraph, integration of the 
above equation uniquely determines Fi{s). 

For simple polygon models such as Ferrers diagrams, the correction-to-scaling functions 
can be written as a polynomial in the scaling function and the scaling variable, as we will 
see below. 
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4 Numerical methods 



We now explain our methods for extracting scaling behaviour of polygon models from given 
series data. The first method involves a numerical analysis of the singular behaviour of 
the area weighted moments by a fit of series enumeration data to the expected asymp- 
totic form. This may then be compared to predictions about the scaling function and 
corrections to scaling. Another independent method consists in the computation of q- 
algebraic approximants. Here, series enumeration data is used to compute an approximate 
g-algebraic functional equation. For polygon models, we believe this to be a good approx- 
imation about the critical point, since the class of g-algebraic functions is rather general 
and since the scaling function seems to be determined by the nature of the singularity of 
the perimeter generating function. The method will then result in approximate differential 
equations for the scaling functions and will provide accurate estimates of its parameters. 
We focus here on algebraic singularities, which is expressed in the requirement that the 
perimeter generating function shall be approximated by algebraic functions [Hj. 



4.1 Moment analysis 

We show how to estimate amplitudes of the area weighted moments from given series 
enumeration data. This in turn provides an estimate of the scaling function and correction- 
to-scaling functions, as follows from paragraph (2.4). 

The coefficients fk^i of the area weighted moments ()2.42|) determine the asymptotic 
growth of the Taylor coefficients of the functions gk{x). This can be seen from the fact 
that asymptotically 



(Xc — 

n^f^^ 7(7 - 1) ^ 7(7-l)(7-2)(37-l) ^ _3^' ^^'^^ 
x« V 2n + 24n2 + ^^"^ > 

The analysis of series data is typically done by a direct fit to the expected asymptotic form 
[inillZI. This provides estimates of amplitudes A^^i 

which can be used to obtain estimates of the coefficients fkj- In particular, we have 

^^,0 = Jp' y (4.3) 

Xc r(7fc,o) 

This expresses the coefficients in the asymptotic expansion of the scaling function with the 
leading asymptotic behaviour of the area moment expansion. The subsequent amplitudes 
determine the functions describing the corrections to scaling. For these amplitudes, the 
contributions from each term in ()2.42|) mix and have to be analyzed carefully. The usual 
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procedure is to assume a special sequence of exponents e(l),e(2), . . . and to analyze the 
rate of convergence of the approximation. If exponents are omitted, then the coefficients 
in the approximation will tend to infinity. If assumed exponents are not present, their 
estimated amplitude will tend to zero. A more detailed exposition of the method may be 
found in [ElITIj. 

4.2 g-algebraic approximants 

The idea of g-algebraic approximants is to fit a g-algebraic functional equation to a given 
function G{x, q) 

M 

P{yo, - ■ ■ ,yN,x,q) = ^^ak,„(a;,g)l/fci---l/fc„ + ao(x,g) = (4.4) 

m=l km 

such that ()4.4|) is exact up to a given degree in the variables. We restrict ai^^{x,q) and 
ao{x,q) to be polynomials in x and q of degrees not exceeding and dg, respectively. 
The coefficients of the polynomials can be found by solving the system of linear equations 
deriving from the expansion of ()4.4|1 in its variables. This process is not unique as there 
are different possibilities to order the set of equations. We demand that the limit g — 1 of 
the approximation leads to algebraic approximants |3] to the perimeter generating function 
G{x, 1). This is achieved as follows. We fix a maximal degree rix- We expand ()4.4|) about 
X = and g = 1 up to the maximal degree rix- This yields a finite number of equations. 
We order this set of equations in increasing powers of 1 — g and, for each power, by 
increasing powers of x. The approximant is obtained by taking, for each power of 1 — g, 
as many equations as possible to obtain a nontrivial solution of the resulting system. This 
gives the g-algebraic approximant of degree dx and dg. In the limit g — > 1, it reduces to the 
algebraic approximant of degree d^ in x. If the class of g-algebraic approximant is correctly 
chosen, the approximate scaling functions will converge to the correct scaling behaviour 
with increasing precision of approximation d^ and dg. 

g-linear approximants have been introduced recently |28j , together with a discussion of 
the relation between this new type of approximants and the existing method of partial dif- 
ferential approximants jHI]. g-linear approxiamnts of first order are suited for the analysis 
of scaling behaviour about poles. They led to the detection of the exact scaling functions 
for exactly solvable g-linear polygon models with a complicated functional equation of 
higher order |28j. We will demonstrate our approach for g-quadratic approximants. 

We have shown above that a square-root singularity leads to a scaling function of 
Riccati type. Since this behaviour is already exhibited by a g-algebraic functional equation 
of degree two, g-quadratic approximants of first order may be employed to approximate 
such a scaling function. Note however that the corrections to scaling may not be correctly 
approximated since contributions from higher degrees of the functional equation enter here. 
As a g-quadratic approximant we choose 

P{yo, yi, X, q) =yl + a{x, q)yoyi - b{x, q)yo + c(x, g). (4.5) 
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From the considerations above, the critical point Xc is given by 



b\x„ 1) = 4c(x„ 1) (1 + a(xe, 1)) . (4.6) 

The coefficients of the differential equation for the scaling function are 

2 _ a'b^ b'b c' 

~ 4(l + a)3 2(l + a)2 + (1 + a)' ^ ' 

^b 

-4/i,o = 



2(l + a)2^ 



where the prime denotes differentiation, and all functions are evaluated at (x, q) = (xc, 1). 
We will see below that this choice of approximant converges well for column-convex poly- 
gons and rooted SAPs, which both have a square-root singularity in their perimeter gen- 
erating function. 



5 Examples 

We first discuss exactly solvable examples, which illustrate the general considerations 
above. These are rectangles, Ferrers diagrams and staircase polygons. Two other exactly 
solvable examples are convex polygons and column-convex polygons, which obey compli- 
cated g-algebraic functional equations. We use these models in order to test our numerical 
tools. We finally analyze the model of self-avoiding polygons numerically. 

5.1 Rectangles 

The perimeter and area generating function of rectangles satisfies the g-linear functional 
equation 

P{yo,yi,x,q) = gx^(l - qx)yi - (1 - qx)yo + qx'^{l + qx) = 0. (5.1) 

It can be derived from the expression given previously |24j for the anisotropic perimeter 
and area generating function by symmetrization or by a graphical argument similar to that 
given previously [21]. The perimeter generating function is rational with a double pole at 
the critical point {x,q) = (1,1). Using an argument analogous to the one given above 
for a simple pole, we assume that the perimeter and area generating function admits an 
expansion about the critical point of the form 

oo 

G(x,g)=e-i^etFfc(se-^), (5.2) 

A:=0 

where s = Xc — x and e = 1 — q. It is seen by arguments analogous to those above that the 
expansion of the g-linear functional equation yields differential equations for the functions 
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Fk{s). The first few functions are given by 



(5.3) 

2) • ' 2 12 12 6 



where Ei(x) = f^e~^^/tdt denotes the exponential integral. The coefficients fn,o in the 
asymptotic expansion of the scaling function are given by fn,o = (— l)""^!- All corrections 
to scaling appear to be linear in Fq{s). Despite the simple form of the functional equation, 
a recurrence for the coefficients is not transparent. We note in passing that the scaling 
function is singular at the origin. This is possibly related to the fact that the condition 
f|2.33|) is violated for rectangles. 

5.2 Ferrers diagrams 

The perimeter and area generating function of Ferrers diagrams satisfies the g-linear func- 
tional equation 

P{yo, yi, X, q) = qx^yi - (1 - qxfy^ + qx^ = 0. (5.4) 

This equation can be derived from the functional equation for the anisotropic model given 
previously |21] by symmetrization or, alternatively, by a graphical argument similar to 
that given there. The perimeter generating function is rational with a simple pole at 
X = 1/2. Thus the considerations of the previous section apply. About the critical point 
(x, q) = (1/2, 1) we assume an expansion of the perimeter and area generating function of 
the form 

00 

G'(x,g)=e-^^eiFfc(se-^). (5.5) 

The scaling function Fq{s) has first been given previously j2H] and has also been extracted 
for a semi-continuous version of the model j2lj. The first few functions Fk{s), as computed 



Fo{s) = Ei{s')e'' 

F,{s) = s^Fo{s)-s-s-' 

^ ' V2 12 2 

' VS 12 3 
1 3 1 1 1 
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from the functional equation, are 



Fo{s) = y|erfc(v^s)e«^' (5.6) 

F2{s) - {—s +—S -As --jFo{s)--s --S + -s 

,16384 „ 1024 7 512 , 124 , 3 \ , 

F3 s = s^ + + s^ s^--s]Fo(s) 

A ^ ^81 3 15 9 A J ^ ^ 

2048 „ 3328 . 256 . 427 . 7 
s - —^s^ - —s^ + —s' + — , 



135 270 135' 

where erfc(x) is the complementary error-function. All corrections to scaling are linear 
in the scaling function, with polynomial coefficients in the scaling variable s. A simple 
recursion for these polynomials is not apparent. 



5.3 Staircase polygons 

The perimeter and area generating function of staircase polygons satisfies the g-quadratic 
functional equation 

PiVo, 2/1, X, q) = yiVo + {2xq - l)yo + x'^q = 0. (5.7) 

It can be derived from the functional equation for the anisotropic perimeter and area gen- 
erating function given previously [22] by symmetrization. Alternatively, it can be derived 
from a concatenation argument similar to that given there. The perimeter generating func- 
tion has a square-root singularity at Xc = 1/4, such that the considerations of the previous 
section apply. About the critical point (x, g) = (1/4,1), we assume an expansion of the 
perimeter and area generating function of the form 

00 

G{x,q) = G^''3\x,q) + ^e^Ffc (se"!) , (5.8) 

fc=0 

where s = Xc~ x and t = 1 — q. The scaling function Fq{s) and the correction-to-scaling 
functions Fk{s)^ where > 0, can be computed recursively from the expansion of the 
functional equation. The scaling function Fq{s) was first computed |22.^ using an explicit 
expression for the perimeter and area generating function. It was recently extracted from 
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the g-functional equation [21]. The first few contributions Fk{s) are 

— -^logAi (2h^ 
16 ds ^ V y 

32 



Fois) 
Fiis) 

F2is) 



Fsis) 
FJs) 



1 d 


— - —sFo{s) + —s^F'is) 
40 15 ^ ' 15 ' 

Ifo{s) - Isf;,{s) 



89 



350 



1936 
1575 



s^Fo(s) + 



525 



2240 



F'(s) + 



(5.9) 



Note that the first correction to scahng vanishes. The above expressions suggest that all 
corrections to scaling may be expressed in terms of derivatives of the scaling function with 
polynomial coefficients. A simple recursion for the polynomials is not transparent. 



5.4 Convex polygons 

The anisotropic perimeter and area generating function G{x,y,q) = C{x) of convex poly- 
gons [ini C] satisfies the following system of g-algebraic equations 



xyq 




xq 
1 



B 



(xq) 



(5.10) 



2xq x^q 
1 + xq xq 
2 1 




where the functions V and W are given by 

xyq 



V = B 



1 — xg' 



W = Bi 



xyq 
1 — xq 



(5.11) 



The function A{x) is the anisotropic perimeter and area generating function for staircase 
polygons. This system of g-algebraic equations can be reduced to a single g-algebraic 
functional equation involving only G{x,y,q): Shifting the parameter x to qx in the above 
system gives five new equations, but only four new variables A{q'^x), B{q'^x), Bi{q'^x), 
Ci(g^x) and C2{q^x). This process is repeated as long as the number of variables remains 
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larger than the number of equations. The variables in this system of equations can then all 
be eliminated, leading to a g-algebraic functional equation for the anisotropic generating 
function G{x,y,q). Since the generating function is symmetric, a g-algebraic generating 
function for the isotropic model G{t,t,q) can in principle be derived from that of the 
anisotropic model by symmetrization. It is however very complicated. 

We restrict ourselves to the perimeter and area generating function of the isotropic 
model in the following. We used the above system to generate series data for the anisotropic 
perimeter and area generating function. We analyzed the isotropic case y = x and observed 
that the area moments go{x) to gji^x) satisfy an algebraic equation of degree four. The 
perimeter generating function jl] is given by 

, - 6a; + llx^ - 4x^) Ax"^ i'cion 

^y""^^)- (l-4x)2 (l-4x)3/2- ^^-^^^ 

It has a singularity at x = 1/4 with a double pole as leading singularity. We conjecture 
that the singular behaviour of all area moments is of the form 

- (l_4a;)4n+2 + (l_43,)3n+3/2 ^^-^"^^ 

/ ^ _ ^^2^+1 (a;) S2n+l(x) 
g2n+l[X) — J- 737rT7 + 



;i_4x)4n+4 (1 _4a;)3n+5/2' 

where the polynomials r„(x) and s„(x) are non-zero at a; = Xc- The leading singular terms 
give rise to an asymptotic expansion of the scaling function 

111 3 

jp (g) ^ ^ „-2 _ ^ -4 , ^ -6 -8 

' 2048 65536 1048576 33554432 

+ s-io s-'^ + s-'^ (5.14) 

268435456 8589934592 137438953472 

4398046511104 ^ ^ ' 

where s = Xc — x. This we recognize as 

"64 



e 

Fo{s) = — — Ei(32s^), (5.15) 



where Ei(s) is the exponential integral. The scaling function for convex polygons is the 
same as for rectangles. This means that convex polygons scale in a first approximation 
like rectangles. The first correction to scaling due to the simple pole is linear in the scaling 
function, as in the rectangle case. We have not found a closed expression for the leading 
corrections to scaling due to the square-root singularity. 
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5.5 Column-convex polygons 

The perimeter and area generating function of column- convex polygons G{x, y, q) satisfies 
the following system of g-quadratic functional equations fS] 

A{x) =(1 + B{qx) — A{qx)){qxy + qxB{x)) + yA{qx) + 2A{x)A{qx) 
B{x) =A{x) + {C{qx) - B{qx))qx{y + B{x)) + {y + 2A{x))B{qx) 
C{x) =B{x) + {C{qx) - B{qx)){qxy + qxC{x)) + 
+ iy + 2B{x))B{qx) + yC{qx), 

where G{x,y,q) = A{x). This system of g-algebraic equations can be reduced to a single 
g-algebraic functional equation involving only G{x, y, q): Shifting the parameter x to qx in 
the above system gives three new equations, but only two new variables B{q'^x), C{q'^x). 
This process is repeated as long as the number of variables remains larger than the number 
of equations. The variables can then all be eliminated, leading to a g-algebraic functional 
equation for the anisotropic generating function G{x,y,q). 

We consider the behaviour of the isotropic perimeter and area generating function 
G(t,q) = G(t,t,q). Since the generating function is not symmetric in x and y, we can- 
not conclude from the above argument that the isotropic generating function G{t, q) is 
g-algebraic as well. The limit g — 1 in (j5.16p yields an algebraic equation of degree 
four for the perimeter generating function Its singularity of smallest modulus is of 
square- root type at tc = 3 — 2-\/2. Assuming an underlying g-algebraic functional equation, 
we can derive the exact scaling function about (t, q) = (tc, 1) by extracting the leading 
singular behaviour of the perimeter generating function and its first area weighted mo- 
ment^. These two functions satisfy algebraic equations of degree four p. The amplitude 
/o = 0.13288910 ... is given as a root of the quartic equation 

4879681t^ + 52908512t^ - 6131584t2 - 161792t + 4096 = 0. (5.17) 

The amplitude /i = —0.0035579552 ... is given as a root of the quartic equation 

565504t^ - 1580288t^ + 216672^^ + 1072t +1 = 0. (5.18) 

Assuming that the scaling function is given by ()3.9|) . we can compare the predictions for the 
leading singular behaviour of the area moments with numerical data from series extrapola- 
tion. More precisely, we compare the amplitudes obtained by numerical analysis with those 
predicted by the scaling function ()3.9p . We therefore generated the series expansion of the 
first ten area weighted moments from the functional equation up to perimeter 40 and used 
the method described above. We assumed corrections to the asymptotic behaviour with 
half-integer exponents, e(m) = m/2, and performed an approximation with ten correction 
terms. The amplitude estimates are shown in the table below. 

^The scaling function for a semicontinuous version of the model has been derived earlier [S]. 
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xlO" 


-8 





516(1) 


xlO' 


"7 


8 


0.3006(1) 


xlO 


-9 


-0.6037(1) 


xlO" 


-9 


-0 


122(1) 


xlO^ 


-7 


9 


-0.7093(1) 


xlO^ 


10 


0.1422(1) 


xlO 


-9 





330(1) 


xlO- 


-8 


10 


0.1877(1) xlO" 


10 


-0.3677(1) 


XlO" 


10 


-0 


997(1) 


xlO- 


-9 



The values /o,i and /i,2 are zero up to numerical accuracy. They corresond to non- 
negative integer exponents, which cannot be extracted from finite size extrapolation, since 
they are regular terms in the corresponding generating function. A comparison with the 
expressions fk^ computed from the assumed exact scaling function coincides with the data 
up to numerical accuracy. We therefore assume that the expression given for the scaling 
function is correct. 

In contrast to the model of staircase polygons, the first correction to scaling is non- 
vanishing. An analysis along the lines of that given in section (3.2) yields the linear 
differential equation 

F[{s) + CoF,is)Fois) = ci + C2sFois) + C3F^'(s) (5.19) 

for the first correction to scaling -Fi(s) with constants c^, where k = 0, . . . ,3. We have in 
particular cq = — (2/i^o)~^- A numerical analysis now reveals that the constant C3 is several 
orders of magnitude smaller than ci,C2. Assuming C3 = 0, we can write solutions of the 
above equation in the form 

F2is) = aos + aiF^(s) (5.20) 

with constants oq, ai. We may impose oq = 0, since any contribution may be absorbed in 
the regular part. The constant oi = —0.0414654242 . . . may be extracted from the defining 
equations of the perimeter generating function and the area moment. It is given as a root 
of the equation 

12032t^ - 67328t^ + 88288t^ + 8432t + 193 = 0. (5.21) 

Using this Ansatz, the values of the consecutive area moments are predicted correctly 
within numerical accuracy. We therefore assume that the above expression is correct. It is 
in principle possible to derive algebraic equations for greater moments of the area to check 
consistency. This requires much higher computational effort, however. Higher corrections 
to scaling may also be analyzed using the above methods. A numerical test is however more 
difficult to perform accurately because the number of undetermined constants increases. 

We also applied the method of g-quadratic approximants in order to obtain estimates 
for the amplitudes /o,o and /i^o- To this end, we computed approximants of the form ()4.5|) 
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to the series data for polynomials of degree one in e and up to degree 9 in t. We get 
/go = 0.1328(1) and /i_o = —0.003557(1), agreeing with the exact values up to numerical 
accuracy. 

5.6 Self-avoiding polygons 

The problem of self-avoiding polygons fSlElEl is an as yet unsolved combinatorial prob- 
lem. The asymptotic form of the finite-size scaling function in the regime g > 1 has been 
determined recently |2E]- It has recently been proved that the anisotropic perimeter gen- 
erating function is not D-finite j2Zj. This suggests that the isotropic perimeter generating 
function has the same property. If so, it cannot then be algebraic. 

Extensive numerical analysis of series data has shown that the perimeter generating 
function has a finite singularity with exponent 3/2 at the critical point Xc = 0.143680628(1) 
|13| ITT)] . This means that the perimeter and area generating function for rooted SAPs 
G^^\x,q) = x-^G{x,q) has a square-root singularity. About the critical point, numerical 
analysis of series data JH] strongly supports an expansion of the form ()2.16p with exponents 
increasing by one (instead of 1/2). We previously analyzed series enumeration data for 
SAPs on the square-lattice, counted by perimeter and area, and extracted the functional 
form of the scaling function j2H]- Here, we explain this analysis in more detail and extend 
it to the first two corrections to scaling. 

We analyzed series enumeration data of the first ten area-weighted moments of the 
half-perimeter generating function for rooted SAPs. The data was produced up to half- 
perimeter 43 using the finite-lattice method and is available on request [23 • We found 
the same exponents = 2/3 and 6' = 1/3 as for staircase polygons, confirming a previous 
investigation [22] • We then assumed corrections to the asymptotic behaviour with half- 
integer exponents, e(m) = m/2, and performed an approximation with ten correction 
terms. The first few values of the amplitudes fk,o, fk,2, fk,3, obtained by an analysis similar 
to that described in the previous section, are shown in the table below. 



k 


fk,0 


fk,2 


fk,3 





-0.9296(1) 


xlO 


-0 


-2.81(1) 


xlO- 


-0 


0.00(1) 


xlO 


-0 
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xlO 
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xlO- 
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xlO 


-0 
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xlO 


-4 
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xlO- 


-3 
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xlO 


-2 


3 


-0.3243(1) 


xlO 


-5 


0.487(1) 


xlO" 
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0.13(1) 


xlO 


-3 
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0.1836(1) 


xlO 
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-0.386(1) 


xlO 


-5 


-0.61(1) 


xlO 


-5 
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-0.1386(1) 


xlO 


-7 


0.374(1) 


xlO- 


-6 
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xlO 


-6 
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0.1301(1) 


xlO 


-8 
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xlO" 


-7 


-0.35(1) 


xlO 


-7 
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-0.1460(1) 


xlO 


-9 


0.568(1) 


xlO" 


-8 


0.37(1) 


xlO 


-8 
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0.1904(1) 


xlO" 


10 


-0.856(1) 


xlO" 


-9 


-0.47(1) 


xlO 


-9 
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-0.2832(1) 


xlO" 
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0.144(1) 


xlO" 
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0.67(1) 


xlO^ 


10 


10 


0.4729(1) 


xlO^ 


12 


-0.269(1) xlO" 


10 


-0.10(1) 


xlO- 


10 



We also determined the coefficients fk,i for the first ten area weighted moments. They 
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are all several orders of magnitude less than the amplitudes fk,o, fk,2, fk,3- Therefore, we 
assume that they are all zero and do not list them in the table. It has been argued by 
universality [7| that /i^o is given by —Xc/ (Stt). The universal amplitude combinations given 
previously |22j may be reconstructed from the values fkfi using ()4.H|1 . 

Though the perimeter generating function is almost certainly not D-finite, it is still 
possible that rooted SAPs are described by a g-algebraic functional equation, leading to an 
algebraic differential equation for the perimeter generating function at some order e" with 
n > 0. Under this assumption, it is still possible to perform a scaling analysis along the 
lines presented above. In fact the argument is exactly the same as that presented above, so 
we conclude that we might expect the differential equation ()3.7|1 for the scaling function. 
To test this assumption, we take the numerical estimates /o,o and fi^ and compare the 
predictions for the amplitudes o for A; > 1 to the numerical estimates in the table 
according to ()2.23p . The prediction coincides with the amplitudes for k > 1 with an error 
less than 10~^. We conclude that the assumed form of the scaling function (|3.9p for rooted 
SAPs is correct. Universality implies that the functional form of the scaling function does 
not depend on the underlying lattice. Therefore, the coefficients c„ of ()2.23|) should be 
obtainable on any lattice from the amplitude combinations 

Cn = fnfif 0,0^ f 1,0 = i~^Tfn,ofofi^ ^'x~^ ' (5.22) 

up to normalization due to the given geometry. This has been successfully demonstrated 
for the triangular lattice PU] . 

It is also possible to apply the method of g-quadratic approximants in order to obtain 
estimates of the constants /o,o and /i^o- We get /o,o = —0.9296(1) and /i_o = —0.005716(1). 
These values coincide with the values obtained by moment analysis within numerical ac- 
curacy given by the available series data. 

We now analyze the corrections to scaling. In order to obtain information about the 
first correction to scaling -Fi(s), we have estimated the coefficients fk,i for the first ten area 
weighted moments. As mentioned above, they are all several orders of magnitude less than 
the terms fkfi and fk,2- We therefore conclude that the first correction is absent, that is 
Fi{s) = 0. 

Under the assumption that rooted SAPs satisfy a g-algebraic functional equation, power 
counting and use of the differential equation for the scaling function results in a differential 
equation for the second correction to scaling F2{s). It is given by 

F^{s) + coF2{s)Fo{s) = cis" + 02^0(5) + c^sF^{s) + c^F^^s) (5.23) 

with constants c^, where k = 0, ... ,4. Remarkably, the solution of this equation can be 
written as a linear combination of derivatives of the scaling function 

F2(s) = ao + ais^ + a2sFo(s) + a^s'^F^is) + a^F^^s) (5.24) 

with constants a^, where /c = 0, . . . ,4. We used this Ansatz for -^2(5) and determined the 
constants from series enumeration data. Note that knowledge of the constants 02, as 
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and a4 predicts the amplitudes of all higher area moments. Again, a numerical analysis 
reveals that the coefficient 04 is several orders of magnitude smaller than 02 and 03. We 
therefore assume that it is equal to zero. We extracted the numbers 02 and 03 by solving 
the Ansatz with series data at order s"^/^ and yielding 02 = 1.01(1) and = 4.01(1). 
The other numbers are then given by ao = —0.01(1) and ai = 0.00(1). The ffist corrections 
to the area moments are then predicted with a relative error less than 10~^. Therefore, 
the above analysis leads to the conjecture 

F2{s) = sFo{s)+4s^F;,{s). (5.25) 

It is difficult to analyze higher corrections to scaling since the numerical accuracy in the 
determination of the amplitudes fk^i decreases with increasing /. 

The scaling function and correction-to-scaling functions for SAPs can be obtained from 
the corresponding functions for rooted SAPs by integration. 

6 Conclusion 

We have analyzed the singular behaviour of g-algebraic functional equations which describe 
exactly solvable polygon models. We considered the special case where the limit g — >• 1 
leads to an algebraic equation. We analyzed the two simplest types of singularity, a square- 
root singularity and a simple pole. It was observed that the scaling function obtained is 
universal - it does not depend on the degree of the functional equation. We gave simple 
exactly solvable examples for each type of singularity. It is obvious that the methods 
introduced may be used to analyze other types of algebraic singularity in detail. For 
example, a finite singularity of type (xc — x)^/" leads to exponents = n/{2n — 1) and 
9 = l/(2n — 1), resulting in a differential equation for the scaling function of the form 
F{s)"' + aF'{s) + bs = 0. This generalizes the result for a square-root singularity and 
may be compared to exact scaling functions predicted recently 0. A related question is 
which polygon models realize these more complicated types of singularity? A natural class 
to look at are models of interacting polymers, whose exactly solvable members are again 
described by g-algebraic functional equations 

Although we discussed only the simplest forms of scaling behaviour, this led to an 
apparently correct prediction for the unsolved model of self-avoiding polygons, where we 
numerically conffimed predictions of the functional form of the scaling function and its 
ffist two corrections to scaling. We know however that anisotropic (rooted) self-avoiding 
polygons are not D-finite j^Zj. This suggests that rooted self-avoiding polygons may have a 
g-algebraic perimeter and area generating function, a conjecture which would be interesting 
to consider more deeply. 

A further step in the analysis of g-algebraic functional equations consists in going be- 
yond the algebraic case, resulting in a differentiably algebraic perimeter generating func- 
tion. Although the analysis in the most general case is difficult, since the singularity 
structure of these types of functions is not fully understood, it is possible to analyze at 
least the case of a linear differential equation in some detail. This case arises from the 
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g — s> 1 limit of g-linear functional equations. This may shed light on the scaling behaviour 
of polygon models such as three-choice polygons ^1] , whose perimeter generating function 
diverges logarithmically. 

From a mathematical point of view, the scaling behaviour of g-algebraic functional 
equations provides examples of uniform asymptotic expansions. There are only few cases 
known where the existence of such an expansion has been proven, for example for staircase 
polygons [22] 5 using an explicit expression for the perimeter and area generating function. 
In the present paper, we extracted the scaling function and corrections to scaling assuming 
the existence of such an expansion and verified the results numerically. It may be possible 
to prove existence directly from the defining g-functional equation, which would answer 
the problem in a much more general framework. 

Self- avoiding polygons are related to a number of other systems via the 0(n)-mapping 
jS2] • This poses the question whether the methods and results of scaling analysis developed 
so far have applications to percolation and magnetic systems, where only a few scaling 
results are available. 
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